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Q Abstract. An enhanced electroweak penguin amplitude due to the presence of unknown new 

I . physics can explain the discrepancies found between theory and experiment in the B — y nK 

decays, in particular in Acp(B — > ty°K ) — Acp(-B° — ► tt + K ), but the current precision of 
the theoretical and experimental results does not allow to draw a firm conclusion. We argue 
that the B a — > <f>p° and B s — s- (f>ir° decays offer an additional tool to investigate this possibility. 
These purely isospin-violating decays are dominated by electroweak penguins and we show that 
in presence of a new physics contribution their branching ratio can be enhanced by about an 
^ ■ order of magnitude, without violating any constraints from other hadronic B decays. This 

makes them very interesting modes for LHCb and future B factories. In [T] we have performed 
J- J^J ■ both a model-independent analysis and a study within realistic New Physics models such as 

' a modified- Z°-penguin scenario, a model with an additional Z' boson and the MSSM. In this 

• , article we summarise the most important results of our study. 

o 



1. Introduction 

The four B — > irK decays have been very useful in testing flavour structure and CP violation 
in the Standard Model (SM) since the late 1990s. These decays have a small branching 
ratio, 0(1O -6 ), and are sensitive to New Physics (NP) contributions. In the past years some 
discrepancies between B — > ttK measurements and SM predictions have occurred, provoking 
speculations on a "B — > ttK puzzle" . To date, the measurements of the branching fractions have 
fluctuated towards the SM predictions, the latter still suffering from large hadronic uncertainties, 
and only the following difference of CP asymmetries shows an unexpected behavior [2]: 

A^cp = A CP (B- -> ir°K-) - A CP (B° -> tt + K~) S = 1.9+IJ % e = P ' (14.8 ± 2.8) %. (1) 

The theory result, derived within the framework of QCD factorisation (QCDF) |3j, and the 
experimental result, taken from [3], differ by 2.5a. This result suggests a violation of the 
isospin symmetry beyond the amount expected in the SM. It can be interpreted as a hint of a 
NP electroweak (EW) penguin amplitude, or alternatively as an enhanced low-energy hadronic 
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effect which the theory is not able to catch. Present data do not allow to distinguish between 
these two possibilities. 

Our aim is to show that the question can be partially addressed exploiting the large variety of 
non-leptonic B decays into two charmless mesons. For this purpose, in [U [5] we have studied the 
purely isospin- violating decays B s — > (f)p° and B s — > 4>7T°, which are dominated by EW penguins, 
as pointed out for the first time in [6] . We have shown that if NP in this sector exists at a level 
where it can explain the AAqp puzzle, it could be clearly visible in these purely isospin-violating 
decays. We present here a summary of the results, referring to [T] for a more detailed discussion. 

2. Isospin violation in B — > ttK and B s — > cf>p°, B s — >■ </>7r° decays 

The B — > irK decays are dominated by the isospin-conserving QCD penguin amplitude. 
Nevertheless, they receive small contributions from the tree and the EW penguin amplitude, 
which are isospin-violating. Combining measurements of the four different decay modes 
B~ — > 7T~K°, B~ — > ir°K~, B° — > ir + K~ and B° — > 7r°K°, it is possible to construct observables 
in which the leading contribution from the QCD penguin drops out, so that they are sensitive 
to isospin violation and therefore to a possible new contribution with the structure of the EW 
penguins. 

The decay amplitudes can be written in terms topological amplitudes which are normalised 
to the dominant QCD penguin contribution. Introducing as well a NP EW penguin amplitude, 
such that 

r E w -+ H3W + r E we~ i5 , rg w rg w + ^ w e~ i5 , r£ w -> r£ w + r^ w e~ iS , (2) 
one obtains 

A^4 C P — — 21m (rc) sin 7 + 21m (r EW + fj| w ) sin 5 . (3) 

fC A) 

Here r-^w, r EW represent respectively the SM colour-allowed, colour-suppressed and 
annihilation electroweak penguin amplitudes, while the fg W are the corresponding NP 
contributions. It is clear from ([3]) that a large A^4cp can be obtained both through enhanced 
low-energy effects in Im (rc), as well as by means of a NP contribution. An Im (rc) large enough 
to solve the A^4cp discrepancy, however, would break QCDF, so that within this framework a 
NP solution is favoured. Although one can construct many other observables which are sensitive 
to isospin violation, the colour-allowed EW penguin contribution and the colour-suppressed tree 
contribution cannot be disentangled within the B — > ttK decays alone since they enter the 
amplitudes exclusively in the combination 

fEW - r c e~ n . (4) 
Our proposal is now to get a more complete picture by considering as well the decays 

B s -)■ (j) P ° and B s -> (f)iT Q . 

These decays are purely isospin-violating and their structure in the SM is very simple. Even 
though we are facing the same type of EW penguin vs colour-suppressed tree amplitude pollution 
as in the B — > ttK decays, various elements suggest that the analysis of these decays may be 
interesting: 

• these decays are dominated by the EW penguin amplitude [6], such that a NP effect coming 
from a new EW penguin amplitude should have more spectacular effects in these decays 
than in the QCD penguin dominated B — > ttK: A new EW amplitude of the same order as 
the SM one can enhance the branching ratios by about an order of magnitude. 
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Figure 1: Enhancement factors of the B s — > tftp°,(f)ir branching ratios with respect to their SM values. 
The black dot represents the SM result while the red striped region shows the theoretical uncertainty in 
the SM. The dark green area is the region allowed by the 2 a constraints from B — > ttK^*\ pK (*\ (ftK^ 
and B s — > cfxft, KK decays; for comparison, the light green area represents the area allowed by constraints 
from isospin-sensitive observables only, considering only B — > nK, TrK^*\pK decays. The solid black line 
represents the la CL of the fit with Scp{B° —> ttK°), while the solid grey line represents the la CL of 
the fit without it. Here the scenario qg ^ is shown. 

• the nature of these B s decays is different from the B — > ttK decays, the first being PV, 
VV decays, respectively, while the latter are PP decays. The non-perturbative low-energy 
QCD dynamics is expected not to be correlated among the two classes of decays, and 
the exact relation cannot be determined within QCDF. Phenomenological comparison with 
other decays such as B — > -kit, B — > irp, B — s> pp suggests a decreasing of non-perturbative 
effects when vector mesons appear in the final state. A NP contribution is instead of high- 
energy origin, and its effects can be reliably studied within perturbation theory, allowing 
for a correlation among the B — > irK and the B s decays of interest. 

3. Model-independent analysis 

We first perform a model-independent analysis. We parameterise NP in EW penguins by adding 
corresponding terms Cy'g NP to the Wilson coefficients of the effective weak Hamiltonian, 

Cif p (m w ) = C™(m w ) g 7 % q% = \ q %\e^. ( 5 ) 

Here Cg° denotes the leading-order SM coefficient, the primed operators are obtained from 
the SM Qy^q by flipping the chiralities of the quark fields. The purpose is to perform a x 2_n t 
in order to determine the NP parameters in such a way that they describe better B — s> ttK 
data. In particular we look for a solution of the AAqp discrepancy. Further hadronic decays 
like B — > pK, ttK* , pK* are used to impose additional constraints at the 2 a level. The result 
of the fit is used to study the decays B s — > 07r°,0/9° and quantify a potential enhancement of 
their branching fractions. Such an analysis, correlating different hadronic decay modes, is only 
possible if hadronic matrix elements are calculated from first principles like in the framework of 
QCDF. A method based on flavour symmetries, as it has been used in most studies of B — s> ttK 
decays so far, could not achieve this. In particular, the decays B s — > </>7r°, (ftp are not related to 
any other decay via SU (3) f so their branching fractions cannot be predicted in this way. 




Figure 2: Enhancement factors of the B s — > <fip°,(j)ir branching ratios with respect to their SM values. 
The meaning of the contours and regions is the same as in fig. [1] Here the scenario q' 7 ^ Is displayed. 
The ler region of the fit is the region to the right of the black (grey) curve. 

As far as the B — > irK decays are concerned, the new contributions to the EW penguin 
amplitudes are given by 

E fBW.ie-** = for -<Z 7 ) [(-0.12)1°^ + (-0.02)1°;°^] + 
t=m ' [0-12^1 + 0.021°;°? i] , 

£ ^~ iSi = i*- & [0.10^ + 0.011°;°^] + 

( q9 - q > 9 ) [0.04l°;°J + (-0.005)1°;°^ i] , 
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E ^w,.^' = (v-q'r) [0.03l°;° 7 + (-0.06)1°;^ z] + 
i=7 ' 9 * 7 '' 9 ' (99 - <Z 9 ) [0.0071°;°°^ + (-0.006)1°;°^] . (6) 

From these expressions one draws e.g. the following conclusions: 

• parity-symmetric models with q 7 = q' 7 and qg = q'g do not contribute to B — s> irK at all. 
Therefore such a scenario cannot solve the A^4cp discrepancy. 

• The contributions f EW 7 (/j and f EW 9 (/) tend to cancel each other. Hence in the scenarios 
with q 7 = qg and q' 7 = q'g only a negligible new colour-allowed EW penguin contribution is 
generated. 

• Whereas Re(f^ w (/) ) features the typical colour-suppression with respect to Re(f EW 9 (/)), 

this pattern is not obeyed by the q^ terms. 

• The annihilation coefficient f^ w 7(/) develops a large imaginary part. In scenarios with 

non- vanishing qip this term gives the dominant contribution to A^cp- 
From eq. ^ we see that the A^lcp discrepancy can be solved either through r^w or through 
f EW . Except for parity-symmetric models, any scenario with at least one of q 7 \ different from 
zero can achieve such a solution. The minimal \q\- value needed to reduce the A^Lcp tension 
below the la level varies among different scenarios. One finds e.g. \qj\ > 0.3 for a model with 



only qi / 0, \qg\ > 0.8 for a model with only qg ^ and = \qg\ > 0.4 for a model with 
qj = qg and the primed coefficients being zero. The fact that in the qj = qg case only a small 
NP contribution is needed, in spite of the absence of tewj demonstrates the importance of the 
annihilation term r^ w . Finally, we like to stress that the solution of the AAcp discrepancy via a 
minimal \q\ - value requires the adjustment of the phase 4> to a certain value. Realistic scenarios 
avoiding such a fine-tuning have larger \q\- values, typically \q\ ~ 1. 

For the qv-only and the gg-only scenarios, the result of the fit to B — > ttK data and the 
consequences for the B s decays of interest are shown in figs. figs. 11121 For details about the fit 
and the observables used we refer to [T] . We find that the B — > irK and related decays set quite 
strong constraints on the parameter space, especially in scenarios where qg ^ or q'g ^ 0. This 
basically rules out the possibility of having \qi\ > 5, i.e. NP corrections cannot be much larger 
than the EW penguins of the SM. The SM point is always excluded at the 2 a level as a direct 
consequence of the AAcp data. According to the sign pattern in eq. (JOJ) , the B — > ttK fits of the 
primed and unprimed scenarios in figs. 11121 are related to each other through rotation by 180°. 
The fit works best in the q^ scenario where the best fit point is given by 

\$\ = 1-9 <pf = -100° (+180°). (7) 

This parameter point yields a full agreement of all the B — > ttK observables with the 
experimental mean values. In the q~' = q^ case a plateau of x 2 = points arises due to 
the large theoretical errors. It turns out that the B — > ttK observables are not very sensitive to 
the g^-only scenarios and so the fit does not work well here. 

From figs. 11121 the enhancement Br SM+NP /Br SM of the B s branching fractions can be read off 
with respect to the different constraint- and fit-regions. In most scenarios an enhancement of an 
order of magnitude or more is possible, especially in those involving qlp ^ 0. The fact that large 
parts of the allowed regions do not overlap with the SM uncertainty regions is encouraging. It 
means that, if such NP is realised in nature, it could be possible to probe it easily. 

Exceptions are B s —> (j)7r for q^p = qjp and B s —> 4>lPl (with the subscript L denoting 
longitudinally polarised vector-mesons) for parity-symmetric models, where no enhancement is 
expected. Furthermore, effects in the q' 9 and the q' 7 = q'g scenarios are limited by the small 
allowed region resulting from the B ttK fit. Largest effects occur as expected in the scenarios 
which are least constrained by B — s> irK, i.e. the single q^ and the parity-symmetric models. 
Especially in these cases a B s — > cf)ir measurement would complement B — > ttK data and, 
while the parity-symmetric models lack the motivation via the A^4cp discrepancy, the q' 7 setting 
resolves it with ease. Moreover, we like to stress that B —¥ ttK data alone cannot distinguish 
among opposite-parity scenarios because such scenarios generate equal results for the B — s> ttK 
observables (for 180°-rotated parameter points). Therefore an analysis of B — >■ irK should for 
example be supported by the analysis of a PV decay, suggesting B s — > (f)ir° as an ideal candidate. 

4. Analysis of specific NP models 

The results of the previous section raise the question: which concrete models for NP can provide 
a large new EW penguin amplitude without being excluded by present data? In this section we 
provide a survey of results obtained analysing a set of well-motivated NP models, referring to 
[1] for more details. The main difference with respect to the model-independent analysis is the 
possibility of adding constraints from other flavour processes beyond the hadronic B modes, e.g. 
the semileptonic decay B — > X s e + e~ , the radiative decay B — > X s j and B s -B s mixing. These 
processes usually yield tight constraints on new flavour structures and it has to be investigated 
if the effects in B — > -kK and B s — > (j)p , <p7r° survive these constraints. 




Figure 3: Enhancement factor for the B s — >■ <fip ,<fiir branching ratios with respect to their SM values 
in the modined-Z°-pcnguin scenario. The green area represents the region allowed by the 2a constraints 
from all the considered hadronic decays, while the area inside the dashed blue line represents the region 
allowed by the 2a constraint from semi-lcptonic decays. The areas inside the dashed orange line represent 
the parameter values for which the modified-Z°-penguin would explain the B s -B s mixing phase 



4-1- The modified Z -penguin scenario 

The simplest class of models with large new contributions to EW penguins comprises models 
with a modified Zsb coupling. Such a FCNC coupling can either be generated by integrating 
out new heavy particles, e.g. in supersymmetric models or fourth-generation models, or it can 
exist at tree-level in more exotic scenarios like models with non-sequential quarks. 

The model is described in terms of two parameters, k s £ r , which parameterise the effective 
coupling of the Z° to left- and right-handed quarks, respectively. A scenario with only k s £ ^ 
scenario shares its most important features with the gg-only setup of the model-independent 
study, while the case with only k r ^ is similar to the (/7-only case. This expectation is 
confirmed by the graphs in fig. [3l We only note that we get a 180° rotation due to the signs of 
5Cg and Cj. The main difference is that we now face additional constraints from semileptonic 
decays and B s -B s mixing. The allowed region for the former is given by the interior of the blue 
dashed curve, the allowed region for the latter by the orange areas outside the zone preferred by 
the B — > ttK fit. We see that the A s anomaly of B s -B s mixing cannot be resolved in a modified 
Z scenario when fulfilling at the same time the semileptonic constraints. However, in models 
with loop-induced modified Z couplings the Z exchange contribution to B s -B s mixing has to 
be regarded as a subleading effect and therefore the corresponding constraint does not apply in 
this case. Focussing on the semileptonic constraints, we find that they are compatible with the 
1 a region of the B — >■ ttK fit for all three cases, but draw the FCNC couplings k s £ r to very 
small values. As a consequence, we expect no significant effects in B s — > </>7r°, cj)p . 

4-2. Models with an additional U{\)' gauge symmetry 

The presence of a heavy Z' boson associated with an additional U(l)' gauge symmetry is a 
well-motivated extension of the SM. This additional symmetry has not been invented to solve a 
particular problem of the SM, but rather occurs as a byproduct in many models like e.g. Grand 
Unified Theories, various models of dynamical symmetry breaking and Little-Higgs models. 
The case of interest, in which the flavour effect is assumed to contribute mainly to the 




Figure 4: Enhancement factors of Br(_B s — > 4>ir°) and Br(B s — > 4>p°) for C£ & ^ 0. The red-hatched ring 
corresponds to the SM uncertainty. The green area is allowed by the 2a constraints from all hadronic 
B decays while the regions inside the blue lines are compatible with the constraint from B s -B s mixing. 
From the biggest to the smallest region they stand for £ = 1/10, £ = 1/25 and £ = 1/100, respectively. 



electroweak penguin, can be described in terms of two parameters, Clw corresponding to the 
properly normalised couplings of the Z'-boson to left- and right-handed quarks, respectively. 
The three scenarios with only ££ 6 ^ 0, only ^ and with ££ 6 = £Jg ^ exactly correspond 
to the (77-only, gg-only and the qj = q' 9 scenarios in our model-independent analysis, except for 
a normalisation factor. The exclusion regions from the 2 a constraints and the confidence levels 
from the fit can be immediately read off from figs. [T] and [21 provided one rescales the axes by an 
appropriate normalisation factor and rotates the pictures by 180° to take into account a minus 
sign in the Wilson coefficients. In this scenario no direct constraints arise from semileptonic 
decays, assuming leptophobic Z' couplings. The relation to B s -B s mixing moreover depends 

on an additional parameter £ = 9u ^' ~]~£r-> which is related to the overall strength of the Z' 

coupling to fermions, and to its mass. From the diagrams in fig. H]we see that the B s -B s mixing 
constraint is in general very tight. It prohibits large effects in B s — > </>7r°, cj>p° for realistic values of 
the parameter £ < 1/25, which would correspond for example to gu(iy ~ 9 and Mz' ~ 400GeV. 
We find that enhancement of a factor ~ 5 is possible in the ££ b and Cj^ scenarios whereas no 
effect can occur in the C'l — Cr- For £ = 1/100 the constraints from B s -B s mixing become so 
strong that no effect in B s (/>7r°,0p° would be detectable. A measurement of a significant 
enhancement would therefore set a lower limit on £, equivalent to an upper limit on the Z' mass. 



4.3. MSSM 

The main contribution to the electroweak penguin in supersymmetry comes from Z penguins 
induced by chargino-squark loops. This scenario is equivalent to a modified Z° penguin model 
with only k s ^ different from zero. Constraints from flavour observables and squark masses, 
however, strongly constrain the parameter space available. In Fig. \5\ we provide a zoom of 
the k s £ scenario with a superposition of points obtained from a scan of the MSSM parameter 
space. The constraints from squark and chargino masses and from Br(i? — > X s j) are taken into 
account. One sees that these points are not able to decrease the AAcp discrepancy below the 
2 a level, nor to create a large enhancement in the B s decays branching ratio. 




Figure 5: Zoom of the upper plots in fig. [3] On top we add the ^SM | values resulting from chargino- 
induccd flavour-violating Z couplings in a parameter scan. 



5. Conclusion 

We have studied the possibility of probing isospin- violating NP in hadronic B decays. Our 
analysis is motivated by discrepancies found in B — > irK decays. In particular, the 2.5 a 
discrepancy found in the observable A^4cp can be interpreted as a sign of NP in the EW penguin 
sector of the theory. Since this discrepancy alone is not able to provide evidence of NP, we have 
proposed here to consider the decays B s — > <f>ir° and B s — > (pp° as well, which are extremely 
sensitive to EW penguins. First, we have demonstrated in a model-independent analysis that 
the solution of the B —¥ irK discrepancy is obtained with an additional NP contribution to the 
EW penguin operators Qj\ Q^g of the same order of magnitude as the leading SM coefficient 
Cg M . Then, we have studied the corresponding enhancement of the B s — > </>7r°, </>/0° branching 
ratios for various scenarios, taking into account constraints from other hadronic B decays. The 
results show that in many cases a large enhancement of about an order of magnitude is possible. 
Exceptions are parity-symmetric models, which have no impact on the VV decay B s — > 4>p°, and 
scenarios with (approximately) equal contributions to and Cg , which cancel in B s — > (j)7r°. 

The last part of the analysis concerns a survey of concrete NP models such as modified Z° 
penguin, a model with an additional U(l)' gauge symmetry and the MSSM. In such models 
additional constraints arise from the semileptonic decays B — > X s e + e~ and B — > K*l + l~ 
and from B s -B s mixing. While the solution of the AAqp discrepancy is in some cases still 
possible, like in the modified-Z°-penguin scenario or in a model with an additional Z' boson, the 
enhancement of the B s — > 4>tt°, 4>p° decays is usually strongly constrained by either the constraint 
from semileptonic decays or from B s -B s mixing. Last, we find that the new contribution to the 
EW penguins in the MSSM is always marginal and cannot solve the AAqp discrepancy, nor 
create a large enhancement of the B s decays considered here. 
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